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From pre- to young Planetary Nebulae: radio continuum 
variability 
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ABSTRACT 

Searching for variability, we have observed a sample of hot post-AGB stars and young 
Planetary Nebulae candidates with the Very Large Array at 4.8, 8.4, and 22.4 GHz. 
The sources had been previously detected in the radio continuum, which is a proof 
that the central stars have started ionising their circumstellar envelopes and an in- 
crease in radio flux with time can be expected as a result of the progression of the 
ionisation front. Such a behaviour has been found in IRAS 18062+2410, whose ra- 
dio modelling has allowed us to determine that its ionised mass has increased from 
10 -4 to 3.3xl0 -4 M Q in 8 years and its envelope has become optically thin at lower 
frequencies. 

Different temporal behaviours have been found for three other sources. 
IRAS 17423-1755 has shown a possibly periodic pattern and an inversion of its ra- 
dio spectral index, as expected from a varying stellar wind. We estimate that the 
radio flux arises from a very compact region around the central star (~10 15 cm) with 
an electron density of 2 x 10 6 cm' 3 . IRAS 22568+6141 and 17516-2525 have decreased 
their radio flux densities of about 10% per year over 4 years. 

While a linear increase of the flux density with time points out to the progression 
of the ionisation front in the envelope, decreases as well as quasi-periodic patterns may 
indicate the presence of unstable stellar winds/jets or thick dusty envelopes absorbing 
ionising photons. 
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1 INTRODUCTION 

The evolution of planetary nebulae remains a challenging 
topic in astrophysics. Planetary Nebulae (PN) have been 
observed over a wide wavelength range, from X-ray to radio 
frequencies. Their complex morphologies and the shaping 
mechanisms that produce them are still a matter of debate. 
Companion stars, jets from central stars, magnetic fields, 
dust tori, and interacting winds are some of the possible 
shaping agents suggested as being responsible for various 
PN morph ologies, and a n overlap of their actions cannot be 
ruled out (jKwoklbOOObl ). 

In general, the current theory of PN evolut ion is based 
on th e Interacting Stellar Wind (ISW) model (Kwok et all 
1 19781 ) and its generalised version (|Kahn fe WestJll985l ). The 
ISW model encompasses a high-velocity stellar wind driven 
by the hot central star that runs into the lower-velocity ex- 
panding circumstellar envelope (CSE) remaining from the 
earlier Asymptotic Giant Branch (AGB) phase. However, 
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this model does not account for the shaping, because it as- 
sumes that an asymmetric distribution of matter is already 
present when the wind interaction occurs. Other models take 
into account the possible role of jets, and have been succ ess- 
fully applied to some nebulae (|Sahai fc Traugeil Il998h . It 
has been suggested that the interaction with a companion 
object, even a massive planet, may provide the necessary 
asymmetry (Soker 2006). Also, large scale magnetic fields, 
influencing or determ ining the shapes, mig ht be sustained 
by a dynamo process (jBlackman et alll200ll ). 

The birth of a PN is defined by the formation of an ioni- 
sation front, which is itself a shaping agent and could heavily 
influence the morphology established in earlier evolutionary 
phases. In this context, important information can be pro- 
vided by observations of very young PN and pre-PN with hot 
central stars, where the physical processes associated with 
PN formation are still occurring. The term "young PN" is 
often used in the literature for nebulae like NGC 7027, which 
hosts a central star with T e g ~ 2 x 10 5 K. When we use this 
term in this work, we actually indicate much cooler objects, 
which show in their optical spectra only recombination and 
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low-excitation emissio n lines (see for exam ple the spectra of 
some of our targets in lSuarez et 

3 (|2006l l). 

To investigate the properties of these rare objects in 
transition from the post-AGB to the PN, we have selected 
a sample of pre-PN and young-PN can didates and s e arche d 
for radio emission fro m ionised shells (|Umana et a] (|2004 ); 
ICerrigone eTal (120081 )1. The targets were selected from stars 
classified in the literature as hot post-AGB candidates, 
showing strong far-IR excess and B spectral type features. 
The detection of radio continuum emission is a proof of the 
presence of free electrons and therefore of ionisation. Multi- 
frequency observations allow for a characterisation of the 
origin of the emission (i.e., different spectral indices point 
to different physical conditions). 

The progression of ionisation in the envelope around 
the star can be observed as an increase in ionised mass with 
time, which implies an increasing radio flux density. A few 
objects have been observed frequently enough to allow for 
such an investigation. For example, CRL 618 has been ob- 
served for many years at cm and mm wavelengths and ap- 
pears to sho w time intervals of increasing and decreasing 
flux density i|Sanchez-Contreras et~alll2004r ). SAO 244567 is 
another transition object that shows a decrea se in radio flux 
density over a few years |Umana e t al 2008) instead of an 
increase, as expected on the ba sis of its last 40 years of 
evolution (jBobrowskv et alll999f ) . Several works monitoring 
this class of objects in the H Q indicate that the onset of 
ionisation may be accompanied b y variability whose origin 
is not clear l|Arkhipova et all I2OO1I ) . while variability due to 
stellar pulsatio ns has been found in post-AGB stars of F-G 
spectral types |Hrivnak et al|[2010l ). 

For a better understanding of this particular phase, we 
have observed again all of the targets detected in our previ- 
ous works looking for radio variability. Six new targets fol- 
lowing the same selection criteria of the original sample were 
also observed in this work, after we detected radio emission 
from them in a previous observing run. 



2 OBSERVATIONS 

We observed our targets with the Very Large Array, op- 
erated by the National Radio Astronomy Observatory 
(NRAO), at 4.8, 8.4, and 22.4 GHz. The observations were 
scheduled dynamically and only 5 sources were observed at 
the highest frequency, because of weather conditions and 
limited time available in the dynamic queue in our RA range. 

The runs at the lower frequencies (4.8 and 8.4) were per- 
formed in June 2009 in the CnB and C configurations, and 
during reconfiguration time. Each target was observed for 
about 10 min, preceded and followed by 2 min on a nearby 
phase calibrator. The absolute flux density scale was defined 
by observing 3C48 or 3C286, depending on which of the two 
calibrators was observable when the run was carried out. 
The same flux calibrators were used in the 22.4 GHz runs, 
carried out in September (DnC array) and December (D ar- 
ray) 2009. At the highest frequency, fast switching between 
the target and the phase calibrator was implemented to al- 
low for a proper phase calibration. The on-source time was 
also about 10 min. 

The data were reduced with the Astronomical Image 
Processing System (AIPS), according to the recommended 



reduction process: the data sets were FILLMed into AIPS 
with average opacity correction and nominal sensitivities 
(dowetght—1) , edited to flag interference or any other bad 
points and CALIBrated along with point weights (do- 
caltb—2). Observations performed in different days were re- 
duced separately and antenna positions were corrected with 
the task VLANT. Since both VLA and EVLA antennas were 
included in the array, a baseline-dependent calibration was 
performed for each calibrator (task BLCAL) and the so- 
lutions were then applied in the usual calibration process, 
to avoid excessive closure errors. Models of both 3C48 and 
3C286 were used at all frequencies to properly set the abso- 
lute flux density scale. Maps were obtained using the task 
IMAGR with natural weights and CLEANed performing a 
few hundreds iterations of the CLEAN algorithm. 

The flux density for each source was estimated by fit- 
ting a Gaussian to the unresolved source (task JMFIT), and 
the rms noise was calculated in an area much larger than the 
synthesised beam (>100 beam), without evident sources in 
it (task IMEAN). One source (IRAS 22568+6141) was re- 
solved in two close peaks. In this case the data were first 
tapered to obtain one unresolved source and then a Gaus- 
sian was fitted, to estimate the total flux density. To better 
analyse the single peaks, at 4.8 and 22.4 GHz the Gaus- 
sian fitting was performed on maps obtained with uniform 
weighting (robust— -5) to improve the angular resolutiona 
and distinguish the peaks. In spite of the uniform weight- 
ing, at 22.4 GHz, the peaks overlap; in this case the task 
JMFIT has been used to simultaneously fit two Gaussians. 
The Gaussians thus obtained are centered on the same coor- 
dinates as those at the other frequencies within about 0".2, 
which confirms that the task correctly identifies the compo- 
nents. 



3 RESULTS 

We summarise in Table [1] the flux densities obtained and 
their relative variation per year. The errors reported in 
the table contain a 3% absolute calibration error, a — 
yjrms 2 + (0.03F V ) 2 . We have c ompared the s e new re- 
sults with the values reported in lUmana et all l|2004l ) and 
ICerrigone et all l|2008h looking for variability. About half of 
the targets exhibit stable radio fluxes and only one is clearly 
variable: IRAS 18062+2410. As "clearly variable", we mean 
that at both 4.8 and 8.4 GHz we observe a relative variation 
larger than three times its error. Variability in at least one 
band is observed in six targets, if this criterion is loosen to 
twice the error. Since the variations are only in one band, it 
is unclear whether this is due to observational errors or an 
intrinsic change in the emission. 

For IRAS 22568+6141, which is resolved into two blobs 
of emission, we list in Table[T]the flux densities of its North- 
ern and Southern lobe as well as the emission from the whole 
nebula. As explained in the previous section, the flux density 
of the whole nebula was obtained by fitting a Gaussian on a 
tapered map, where the source is not resolved. The sum of 
the flux densities of the two blobs does not account for the 
emission of the entire nebula. An inspection of the residual 
maps (the maps after subtraction of the fitting Gaussians) 
indicates that about 4.3 mjy at 4.8 GHz and 0.5 mjy at 8.4 
GHz remain unmatched. This means that - besides the two 
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Table 1. Radio flux densities obtained at the VLA and their relative variations per year. Targets that do not show variations are listed 
first, then targets that show variations at least in one band. For IRAS 22568+6141, the flux densities of its Northern and Southern lobes 
are also listed. 



IRAS ID 


F4.8GHZ 


AF 4i8 GHz 


F8.4GHz 


AF 8 .4GHz 


F22.4 GHz 




mJy 


% yr 


mJy 


7o yr 


mJy 


01005+7910 


0.29+0.05 




0.17+0.04 


-4.3+4.9 




17381-1616 


1.64+0.08 


0.3+1.2 


1.62+0.07 


2.3+1.4 




19336-0400 


11.1+0.3 


1.1+1.0 


9.5+0.3 


-0.3+1.0 




19590-1249 


3.3+0.1 


1.1+1.3 


3.1+0.1 


1.6+0.9 


- 


20462+3416 


0.62+0.05 


0.7+2.0 


0.43+0.04 


0.4+2.6 




21546+4721 


1.56+0.06 




1.56+0.06 


0.6+1.4 


1.32+0.07 


22023+5249 


2.60+0.09 




2.28+0.08 


2.3+1.2 


2.1+0.1 


22495+5134 


9.6+0.3 




8.6+0.3 


0.8+1.2 


8.1+0.3 


06556+1623 


0.45+0.05 


-4.1+1.5 


0.48+0.04 


-2.1+1.5 


0.11 a 


17423-1755 


0.63+0.06 




0.50+0.05 


15.4+4.8 




17460-3114 


1.4+0.1 


-1.7+1.4 


1.08+0.07 


-2.7+1.3 




17516-2525 


0.4+0.1 6 




0.32+0.06 


-11.6+2.8 




18062+2410 


1.72+0.08 


4.9+1.5 


2.67+0.09 


13.8+2.2 




18442-1144 


21.7+0.7 


3.7+1.2 


20.7+0.6 


1.3+1.0 




22568+6141 


23.8+0.7 




21.7+0.7 


-8.1+0.7 


16.0+0.5 


22568+6141 N 


10.0+0.3 




11.2+0.3 




8.5+0.2 


22568+6141 S 


9.8+0.3 




9.3+0.3 




7.3+0.3 



a The source was not detected and the rms of the map is listed. 
'Marginal detection. 



peaks - weak extended emission is present and is recovered 
by tapering. 

We performed for the first time multi-frequency ob- 
servations of IRAS 01005+7910, 17516-2525, 21546+4721, 
22023+5249, 22495+5134, and 22568+6141. 

We fitted the data points to derive the spectral indices, 
which turned out to match very well with what expected for 
optically thin shells (a ~ -0.1), as can be seen in Table [2] 
Since the emission is optically thin, we can also derive the 
emission measure for each of these sources as 



< EM >= 



f n EMdQ _ 5.3 x W 5 F 8 aghz 



e 2 



where Fsaghz is the flux density at 8.4 GHz in mJy, an d 
8 the angular radius in arcsec (|Terzian fc Dickevl 1 19731 ). 
The only target that is resolved at our angular resolution 
is IRAS 22568+6141 (§ 0, whose northern and southern 
blobs of emission are listed separately in Tabled Gaussian 
fitting of the two peaks allows us to determine their posi- 
tions as RA=22:58:51.44 DEC=61:57:45.2 (North peak) and 
RA=22:58:51.68 DEC=61:57:42.8 (South peak). The peaks 
are about 3" apart, with PA~146°. 

To calculate their emission measures, we need to esti- 
mate the angular sizes of our sources. In Tabled we list the 
major and minor axes of the convolution beam. For those 
sources that have be en resolved wit h high - angular resolu- 
tion observations in ICerrigone et~al (|200Sl) (IRAS 18442- 
1144, 19336-0400, 19590-1249, and 22023+5249), the geo- 
metric mean of the axes is an approximate estimation of 
their diameters, therefore we take into account this value in 
our calculations. 



Table 2. Spectral indices from data fitting (with statistical er- 
rors), emission measures (lower limit) and beam sizes of the 8.4 
GHz data. The list is divided in sources previously detected (up- 
per list) and new detections (lower list). 



Target 


Spec. Index 


Emission measure 
10 s cm -6 pc 


Beam 
arcsec 2 


06556+1623 


0.10+0.47 


2.6+0.2 


3.1x1.3 


17381-1616 


-0.02+0.22 


2.6+0.1 


4.7x2.8 


17423-1755 


-0.40+0.47 


0.78+0.07 


4.9x2.8 


17460-3114 


-0.51+0.35 


1.30+0.08 


6.9x2.6 


18062+2410 


0.78+0.20 


_a 


3.0x2.7 


18442-1144 


-0.09+0.15 


38.2+0.1 


4.2x2.8 


19336-0400 


-0.27+0.15 


20.4+0.1 


3.6x2.7 


19590-1249 


-0.13+0.17 


9.4+0.1 


3.0x2.3 


20462+3416 


-0.67+0.46 


2.0+0.2 


2.6x1.7 


01005+7910 


-0.93+0.98 


0.8+0.2 


2.8x1.7 


17516-2525 


-0.38+1.29 


0.5+0.1 


5.7x2.2 


21546+4721 


-0.11+0.09 


7.8+0.2 


2.7x1.6 


22023+5249 


-0.14+0.08 


8.6+0.2 


2.6x2.1 


22495+5134 


-0.11+0.06 


11.7+0.2 


5.3x2.9 


22568+6141 


-0.26+0.05 






22568+6141 N 


-0.14+0.02 


148+4 


2.0x0.8 


22568+6141 S 


-0.20+0.04 


123+4 


2.0x0.8 



Not optically thin: see modelling in § \E\ 



4 IRAS 17423-1755 

IRAS 17423-1755 (Hen3-1475) is a well known point- 
symmetric pre-PN at a distance of about 5.8 kpc. It shows 
OH maser emission, lines from ionised elements, and re- 
combination lines due to shocks propagating in its CSE 
llRiera et a 3 120031 ). ISanchez Contreras fc Sahal J2001) find 
that two wind components are present in the vicinity of the 
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Figure 1. Continuum spectra for targets that had not been 
observations. 
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observed before at radio wavelengths or that lacked multi-frequency 



star (~ 0".7): a fast and an ultra-fast wind (150-1200 and 
2300 km s _1 respectively), both with kinematical ages of 
tens of years. Decreasing velocities from 1000 to 150 km s _1 
are found in the knots in its bipolar outflow, traced by op- 
tical and near-IR imaging over 17" (jRiera et a ll l2003h . Mil- 
limeter and cm observations trace much smaller structures, 
likely associated to a circumstellar disk (millimeter obser- 
vations ) and a compact io nised region (centimeter observa- 
tions) ijHuggins et alll2004h . 

Several models have been applied to this source to ac- 
count for its kinematical and morphological properties. A 



possible explanation is that its pecularities are the result 
of an ejection variability of the central source (|Riera et all 
l2003h . I Velazquez et all (2004) have modelled this nebula with 
a precessing jet of periodically variable ejection velocity su- 
perimposed to a linear increase, which propagates into an 
Interstellar Medium (ISM) previously perturbed by an AGB 
wind. The velocity variability in their model has a period of 
120 yr and a half-amplitude of 150 km s _1 . The jet has there- 
fore an ejection velocity of the form Vj et = vo+Visin(2ir(t — 
to)/r) + at, with vq = 400 km s , v\ =150 km s , 
to = — 640yr, t = 120 yr, and a —1 km s _1 yr" 1 . 
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iLee fc Sahail (2003) present a general model to shape 
pre-PN with Collimated Fast Winds (CFW) interacting with 
a spherical AGB wind and they compare their results with 
observations of CRL 618. They introduce periodic variations 
of the density and velocity of the fast wind in such a way 
that the mass-loss rate (M = 4-7rr 2 pv) is kept constant: the 
density varies as p = - 1 ^ P J^ and the velocity as v = 

Vf(l + A sin — ). In particular, to compare their model to 
CRL 618 they consider a period of 22 yr and a velocity half- 
amplitude of 150 km s . 

Radio continuum data since 1991 are available in the 
VLA archive. We retrieved the data, reduced them follow- 
ing the same recipe as for the rest of our data, and found a 
flux density of 0.36±0.04 mjy in 1991 and 0.56±0.04 mjy in 
1993, both at 8.4 GHz. As shown in Figure [21 the combina- 
tion of our 2001 and 2009 observations with the archive data 
clearly shows variability of the radio flux density. Unfortu- 
nately, our temporal sampling is too coarse to conclusively 
assess whether the variability is regular or not. Nevertheless, 
we show that the data are well fitted by periodic functiona l 
for ms similar to those i ntroduced by I Velazquez et a] (|2004l ) 
and lLee fc Sahail (|2003h : 



IRAS 17423-1755 



S v = A sin j^t + A 2 



l + Ax sin - 



(dotted line) 
(dashed line) 



We are fitting four points with functions having three 
degrees of freedom, therefore the procedure should not be 
over-interpreted. Both fits give periods around 19 yr (18.6 
and 18.7 yr, for the dashed and dotted curves, respectively), 
comparable to the periodicity in the model of CRL 618 by 
I Lee fc Sahail (|2003h , but different fro m the period introduce d 
in the model of IRAS 17423-1755 bv lVelazquez et all (|2004l ). 

Our 2009 observations indicate that the emission is op- 
tically thin, with a spectral index of -0.4 (TablelH). Nevert he- 
less, the observations reported in lCerrigone et al (|2008l ) at 
8.4 and 22.4 GHz indicate optically thick emission between 
these frequencies with a spectral index around 1. In the cited 
paper, it was assumed that no variation had occurred in the 
source between 2001 (epoch of the 8.4 GHz data) and 2003 
(epoch of the 22.4 GHz data). Figure [5] shows that this as- 
sumption is not correct and that in 2003 the flux density at 
8.4 GHz was smaller than in 2001, if the variability is reg- 
ular. This does not affect the conclusion that the spectral 
index was positive in 2003 and likely around 1, considering 
errors. 

If the radio emission is due to the stellar wind, the ra- 
dius R of the emitting region depen ds on the frequency as 
R(j/)~ v~ - 7 (|Panagia fc Fellilll975l ). in the optically thick 
regime, as is the case in the 2003 observations. Since any 
variation in the wind will propagate outwards and the emis- 
sion at higher frequency arises from inner regions, a light 
curve at a lower frequency will be delayed compared to one 
at a higher frequency. Therefore, the value of the spectral 
index will also vary with time, being positive when the emis- 
sion is decreasing and negative when it is increasing. The 
measurement of a positive spectral index in 2003 and a neg- 
ative one in 2009 strengthens then the interpretation of the 
observed variability as due to a n unstable stellar wind. As 
shown in iPanagia fc Feul (| 19751 ). if the radio flux is due to 
a stellar wind, we can calculate its mass loss rate as 
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Figure 2. VLA observations at 8.4 GHz of IRAS 17423-1755 over 
16 years. The point in 1991 is from the VLA project AK247 (PI: 
Knapp), while that from 1993 is from AZ63 (PI: Zijlstra). 



M = 0.32 X 10" 



( F^V'* ( D_\ Z/2 



10 3 km s 1 \mjy J \kpc J 
and the size of the emitting region as 

M 



-2/3 



Meyr" 1 



R ^ = 6 ' 23 X 1014 ( 10^) V lO- 5 M yr- 

/ V<x> 

' V10 3 km s" 1 



2/3 



where cosmic abundances, full ionisation, and a 10 K tem- 
perature of the electron gas have been assumed. If we link 
the radio flux to the fast wind component, we can take 
100 km s _1 as its terminal velocity and calculate from the 
radio flux density in 2001 (when the spectral index was likely 
positive) M ~ 1.6 X 10~ 6 M Q yr -1 and R 8 .4GHz ~ 10 15 cm 
(~ 0".01), assuming a distance of 5.8 kpc. 

The mass loss rate found is quite high for a post- AGB 
star, where values smaller than 10 -7 Mq yr -1 are expected. 
Nevertheless, large rates (10 -5 -10 -4 Mq yr -1 ) have been es- 
timated in the prototypical pre-PN CRL 618, which has un- 
dergone at least two distinct episode s of mass loss in the form 
of a s low wind, in the last 2500 yr l|Sanchez-Contreras et all 
l2004h . Also, en hanced mass l oss in IRAS 17423-1755 has 
been found by iHuggins et all l|2004l ). who have estimated 
that its molecular envelope (~0.6 Mq) was ejected in less 
than 1500 yr at a mass loss rate > 10 -4 Mq yr -1 . The esti- 
mation of the linear size of the emitting region allows us to 
improve the calculation of the emission measure and convert 
it into an electron density of 2xl0 6 cm" 3 . 

If we assume that, like in CRL 618, the radio flux arises 
from a slow stellar wind (15 km s _1 ), its mass loss rate 
becomes 2.4 x 10 -7 Mq yr -1 , with Rs.4 approximately un- 
changed. 

The lack of multi-epoch observations covering the whole 
cm range does not allow for conclusiveness. 



5 IRAS 17516-2525 

This source was first studied by IVan der Veen et a 3 d!989h . 
who had selected it as a post-AGB candidate because of its 
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IRAS colours. Hydroxil maser lines with At) ~ 40 km s" 1 
we re detected with the Pa r kes a nd N anqay radiotelescopes 
bv Ite Lintel Hekkert et all l)l99lf l and ISzvmczak fc Gerard! 
l|2004 ) respectively; the l atter detected both 1612 a nd 1667 
MHz lines. Nevertheless, IVan der Veen et all (|l989l ) did not 
detect the OH 1612 MHz line with the VLA in B array at the 
coordinates of the source, but ~7' north of it. It is likely then 
that the detections performed with the single-dish telescopes 
are to be associated to this off-set maser, although further 
interferometric observations would be desirable to rule out 
an i ntrinsic OH variability. 

ISanchez-Contreras et all l|2004l ) performed optical spec- 
troscopy of IRAS 17516-2525 and found that its spectrum 
is dominated by nebular emission lines. They also suggest 
the presence of a hot central star, based on the similarity 
between its spectrum and that of IRAS 19520+2759, which 
they classify as Q9. In their an alysis of the optical spectrum, 
ISanchez-Contreras et all (2004) find that the Ha has a P-Cyg 
profile, indicative of on-going mass loss, although Pa lines do 
not display the same profile. Furthermore, they detect the 
[Call] doublet at 7291.47 and 7323.89 A, which leads them to 
conclude that moderate- velocity (40<v<100 km s" 1 ) shocks 
must propagate in the envelope around this star and distruct 
dust grains, enhancing the amount of Ca in the gas phase, 
which would otherwise be depleted onto grains. 

We detected IRAS 17516-2525 for the first time in 2005 
and find a strong decrease of its flux density over 4 years 
(11% per year), the source having almost halved its emission 
at 8.4 GHz. The spectral index between 4.8 and 8.4 GHz 
indicates an optically thin spectrum, although the detection 
at the lower frequency is only marginal. In the context of the 
observations available in the literature, the radio continuum 
variability may be due to instabilities of the stellar wind (i.e., 
shocks) propagating in the circumstellar envelope. More and 
more sensitive observations are needed to inspect this target 
and determine its radio continuum spectrum, to possibly 
distinguish between an ionised shell and a stellar wind, and 
to look for a variability pattern. Within the uncertainty due 
to the larger error at the lower frequency, our estimate of 
a fiat radio spectral index strengthens the hypothesis that 
this is a young PN. 



6 IRAS 18062+2410 

This source has been monitored at optical wavelengths over 
a wide period of time and shows photometric variability on 
both short (hours) and long (years) time scales. In particu- 
lar, the equivalent widths of H a and Hg increase and its op- 
tical continuum decreases (|Arkhipova et alll2007h . While the 
short-time variability has been explained in terms of stellar 
pulsations and an unstable stellar wind, the long-timescale 
variations have been attributed to an increase of the degree 
of ionisation of the envelope, due to a stellar temperature in- 
creas ing at a rate of about 200 K per year |Arkhipova et all 
l2007h - 

Previous VLA observations of this target perform ed in 
2003 and 2005 are reported in ICerrigone et all IMPOST) . We 
have modelled the radio continuum emssion assuming that 
the central star is surrounded by a shell of ionised gas with 
density radially decreasing as r~ 2 . We have calculated a 3D 
distribution of mass in a spherical shell with outer radius 
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Figure 3. Radio continuum spectra of IRAS 18062+2410 ob- 
tained at different epochs. The correspondent model curves are 
shown in the same colour as the data points. 

Rout, inner radius R;„ = jyRout (with < r\ < 1), and den- 
sity at the inner radius pi n . The calculation of the average 
density in the shell links this parameter to the density at the 
inner radius: < p >~ rj 2 (l — rj)pi n . The electron temperature 
was set to 10 4 K. D etails of the radio model can be found in 
lUmana et al l|200Sl ). We fit the 2003 data with Ri„=0".0354, 
Rout=0".06, pin— 5.2 10 5 cm" 3 , and d=6.4 kpc, which gives 
us an estimated ionised mass about 1.48 10" 4 M . 

The nebula has increased its flux density in 2005 and 
2009, as shown in Figure^ To estimate the expansion veloc- 
ity of the ionisation front, we assume that the CSE expands 
at a typical velocity of 15 km s _1 (therefore its inner radius 
will be 0".0364 in 2005). A fit to the 2005 data enables us 
to calculate that the velocity of the ionisation front is of the 
order of 100-200 km s _1 , the typical range for fast winds in 
post-AGB stars. 

Having the above mentioned values as first guesses, we 
now try to fit all of the 2001-2009 data keeping the CSE 
expansion velocity set at 15 km s _1 . We first set Ri n o, 
Routo (the values in 2001), and p in o (the density at Ri n o, 
which does not vary with the epoch) and v ou t (the veloc- 
ity of the ionisation front within the shell). Then we cal- 
culate R in =R;„o + v in At and R «t=Ro«to + (vin + v out )At 
for the other epochs, At being the time between the first 
observation and the specific epoch; finally, we calculate 
the density at Ki n as pi n =Pm0 ( ) • A satisfactory fit 
can be obtained if v out — 120 km s _1 (kept constant at all 
epochs) and in 2001 we have R;„=0".0380, R olit =0".0522, 
and pi„=5.38 10 5 cm -3 , as shown in Figure[3] 

The models allow us to estimate that the critical fre- 
quency has decreased from around 8.2 GHz in 2001 to 7.6 
GHz in 2009, pointing to an optical thinning of the enve- 
lope, with the opacity at 8.4 GHz ranging between 0.96 in 
2001 and 0.84 in 2009. The ionised masses are 1.0 x 10" 4 , 
1.5 x 10" 4 , 2.2 x 10~ 4 , and 3.3 x 10" 4 M 0) respectively in 
2001, 2003, 2005, and 2009. 

Although a linear increase of the flux density produces 
a good fit to the data, it cannot rule out the possibility of a 
cyclic variability with a period much longer than 6 years. 

Figure 31 shows the variability that we find in this ob - 
ject, which confirms what reported in lCerrigone et al (2008). 
We have overplotted the 8.4 GHz flux density calculated 
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Figur e 4. Radio fl ux density at 8.4 GHz and flux 
l| Arkhipova et all 120071 ) at several epochs for IRAS 18062+2410. 
The dashed curve is a fit to the points before 2002. 

from the fluxes in I Arkhipova et al (2007), which we have 
estimated from the equivalent width s and absolute va l ues in 
1996 and 2006 given in that work. lArkhipova etU \200i ) 
estimate that the Hp flux was 2.7+0.3 10 -13 erg cm -2 s _1 
in 1995-1996 and 6.2+0.4 10" 13 erg cm -2 s _1 in 2006 and 
list its equivalent width over 10 years. To convert the equiv- 
alent width into flux, we have calculated the average width 
over all of the 1995-1996 observations (one single high point 
was neglected) and over the 2006 data, then scaled the av- 
erage values to the respective fluxes given in the paper, and 
finally averaged these two scaling factors together to ob- 
tain one mean equivalent- width/fiux conversion factor. After 
converting into flux units each equivalent width, these were 
converte d into rad i o flux densities at 8.4 GHz following Eq. 
IV-26 in lPottaschl j 19841 ) 

S v = 2.51 x 10 7 T e ' 53 v-° A Y F H[i 

where S„ is the flux density in Jy at the frequency v, T e 
is the electron temperature (set to 10 4 K), v is in GHz, Y 
accounts for contributions from He ions and is set to 1.1, 
and Fh /} is the flux in the H,g in units of erg cm -2 s~ . 

The mismatch between the two data sets in Figure [4] is 
likely due to an optical-depth effect. As we have shown in 
our modelling, the nebula is undergoing optical thinning of 
its envelope. When the radio flux density is optically thick, 
its observed value will be smaller than that calculated from 
the H,g flux; on the contrary, when it is optically thin, it will 
be larger than the value calculated from the H^, because 
of absorption at optical wavelengths. Figure [3] seems to be 
showing the passage from one regime to the other. 

The H^ data seem to show a periodicity before 2002. 
After excluding three clear outsider points, we have calcu- 
lated a fit to the data before 2002 with a functional form 
of type Fua = Aosin^nt / A\) + Ait + A3, which includes 
a periodic term and a linearly increasing one. This fit gives 
a period of ~4.1 yr and a starting time ~6.3 yr before the 
first observation (A =0.47, A x = 1500, A 2 = 8.7172 10~ 4 , 
A3 = -2133.622). 

The number of Lyman photons necessary to sustain the 
ionised region can be calculated as 

Q = VTleCtB 



where V is the volume of the ionised region, n e the aver- 
age electron density, and Qs the recombination coefficient 
summed over all the excited st ates (2.59 x 10 -13 cm 3 s _1 , 
at T e = 10 4 K) jKwokl bOOOal ). Taking advantage of the 
parameters derived from our modelling, we can estimate 
that in 2001 Q = 2.2 x 10 46 photons s" 1 , while in 2009 
Q — 1 .0 x 10 47 photo ns s _1 , which by interpolating the val- 
ues m iPanagial (|l973l ) give temperatures of the central star 
of ~23400 and ~25900 K, in 2001 and 2009 respectively. 



7 IRAS 22568+6141 

IRAS 22568+6141 was discovered by iGarcfa-Lario et all 
jl99ll ). who classified it as a possible young PN and es- 
timated a distance of 6 kpc, electron density about 2 x 
10 4 cm -3 and T c g ~24000 K. For this source, a compari- 
son with previous data is possible only at one frequency. We 
first detected this nebula with the VLA (C array) in July 
2005 at 8.4 GHz and found a flux density of 32.70+0.98 mjy 
(Umana et al, in preparation). The source was also included 
in a survey carried out with the Torun 32-m Radio Telescope 
betwe en December 2005 and May 2007 by IPazderska et all 
(2009). They measured 31.4+4.1 mjy at 30 GHz, which is 
in agreement with our measurement at 8.4 GHz, because we 
find that the emission is optically thin. 

Our observations indicate a decrease in the emission 
around 34% between 2005 and 2009. Since the source is 
known to be extended and we observed it with an in- 
terferometer in different array configurations, missing flux 
could be the reason of the apparent variability. Interfer- 
ometers in fact filter out emission from structures beyond 
a size - the Largest Angular Scale (LAS) - that depends 
on the shortest baseline and on the tw-coverage. Our low- 
frequency observations were carried out during reconfigu- 
ration time and the minimum baselines turned out to be 
about 1.3 and 2.1 kA at 4.8 and 8.4 GHz respectively. These 
correspond respectively to a LAS around 40" and 23" in 
snapshot mode. The high-frequency run was performed in 
D array, with LAS~33". It is therefore unlikely that we are 
missing any flux from extended structures, since this source 
is known to have a size around 8" fro m observations in other 
wavelength ranges (|Sahai et alll2007l ). We can conclude that 
the fading of the nebula is real. 

Figure [5] shows the double-peak morphology found in 
the target, which matches what observed at optical and 
near-IR wavelengths. The southern peak is not resolved in 
any of our maps, while the northern one is partly resolved 
with a size around 2.0" x 0.8" (PA~ 50°). The emission 
measures are 12.3 x 10 6 and 14.8 x 10 6 cm -6 pc, for the 
Southern and Northern lobe respectively, where the same 
size has been used for both blobs. 

If we assume a distance to the nebula of 6 kpc, from 
the emission measure we find densities of 2.0 x 10 4 cm~ 3 
(North peak) and 1.8 x 10 4 cm -3 (South peak). The corre- 
sponding decay times for electron recombination are tdecay = 
(J'/cm -^) yr |Kwokll2000al ). then about 3.8 and 4.2 years for 
the N and S blobs respectively. For the northern peak, which 
is partly resolved in our observations, we also calculate the 
ultraviolet luminosity necessary to produce the observed ra- 
dio flux density as done for IRAS 18062+2410. In this case, 
we calculate the luminosity necessary if the star were at 
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Figure 5. Bipolar morphology observed at 8.4 GHz in IRAS 
22568+6141. The contours are 60x(-3, 3, 9, 27, 81, 100, 
130) /^Jy/beam and the grey scale is in mjy/beam. The con- 
volution beam is shown in the bottom left corner. 



the centre of the blob, assuming the emitting volume is a 
sphere whose radius equals the geometric mean of the semi- 
axes from the Gaussian fit. We obtain Q ~ 4.7 x 10 46 pho- 
tons s _1 . Because of the geometry of the source, this is to 
be regarded as a lower limit to the stellar UV luminosity, 
since the photons ionising the blob are only a fraction of all 
those emitted by the central star, and many may also be ab- 
so rbed by c i rcums tellar dust. By interpolation of the values 
m IPanarial (|l973h for class III, this gives us a lower limit 
of 24600 K for the temperature of the central star, slightl y 
larger than what estimated by iGarria-Lario et all (|l99lf ). 
This lower limit points out a possible rapid increase of the 
temperature of the central star. 

Our observations confirm that this object has started 
its PN evolutionary phase. 



8 CONCLUSIONS 

The transition from post-AGB stars to PN is a crucial point 
of the evolution of intermediate-mass stars. The develop- 
ment of an ionisation front determines major changes in the 
circumstellar environment and may contribute to the shap- 
ing of these intriguing sources. Winds and jets are among 
the agents responsible for the early shaping of the nebulae 
during the early post-AGB phase, if not even before. 

The ionisation of the envelope gives us the possibility to 
study the action of these agents around the star. Observa- 
tions at optical wavelengths have pointed out that hot post- 
AGB stars show photometric variab ility, which in some in - 
stances has a quasi-periodic pattern (|Arkhipova et alll2007l ). 

Within our sample of pre- and young PN, we have 
found that 4 of our targets exhibit radio-continuum vari- 
ability. In IRAS 18062+2410, the variability consists in a 
steadily increasing flux density, while in IRAS 17423-1755 
it may follow a periodic pattern. The remaining two vari- 
able sources, namely IRAS 17516-2525 and 22568+6141, 
have both decreased their flux densities of about 10% per 



year over a 4-year time. The cause of such a decrease is 
not clear. They may follow a variability pattern similar to 
IRAS 17423-1755, but in this source the variability occurs 
in a very compact region around the central star, while in 
IRAS 22568+6141 it involves two regions likely ~1.5" away 
from the star (1.3 x 10 17 cm at 6 kpc). 

In near-IR images, the central star in IRAS 22568+6141 
appears completely obscured, which implies the presence of 
a thick circumstellar region of dust and molecular gas. Even 
at cm wavelengths, we do not detect any central emission, 
although this may be an angular resolution issue. We can 
speculate that this obscuring region is in fact a circumstel- 
lar disk and the two radio blobs are due to the action of a 
collimated wind onto the circumstellar nebula. In this sce- 
nario, if the wind slows down, fewer photons per second will 
reach the CSE. Because of the high density, the decay time 
is short and then the radio flux will decrease, as a result of 
recombinations. In an alternative scenario, we can imagine 
that the central star is rapidly increasing its temperature, 
producing more and more photons at shorter wavelengths. In 
the blackbody approximation, the radiation from the central 
star peaks around 1207 A, if T c ff =24000 K. Down to about 
600 A for silicates and 800 A for amorphous carbon, the 
dust absorption coefficient increases with decreasing wave- 
length, therefore the hardening of the radiation from the 
central star would suffer stronger dust absorption determin- 
ing a decrease of the ionization degree and of the observed 
radio flux density, until the radiation manages to destroy the 
dust grains. A decrease in radio continuum emission over a 
few years might then indicate a fast increase in stellar tem- 
perature. 

Variability may be present in several other targets in 
the sample, although the data do not allow for any con- 
clusive statement. Similar results have been obtained for 
other stars transiting betwe en the late post-AGB and the 
PN phases, such a s CRL 618 JSanchez-Contreras et a]||2004T ) 
and SAO 244567 jUmana et aJl200Sh . 

While a linear increase in time of the flux density points 
out to the progression of the ionisation front in the envelope, 
quasi-periodic patterns may indicate the presence of jets or 
instabilities of the stellar wind (i.e., shocks). Several models 
invoke instabilities and episodic ejectio ns, or quasi-periodic 
variab ility. The MHD simulations in iGarci'a-Segura et al.l 
( 1997) show that the collimated outflow can be subject to 
instabilities and form blobs (i.e., the knots observed in pre- 
PN), or an episodic jet can be produced by periodic varia- 
tions of the magnetic field. Such instabilities typically have 
much longer time-scales than what observed in our targets 
and cannot be linked directly to our (relatively) short-term 
variations. 
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